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Frame dragging (Lense-Thirring effect) is generally associated with rotating astrophysical objects. 
However, it can also be generated by electromagnetic fields if electric and magnetic fields are 
simultaneously present. In most models of astrophysical objects, macroscopic charge neutrality is 
assumed and the entire electromagnetic field is characterized in terms of a magnetic dipole component. 
Hence, the purely electromagnetic contribution to the frame dragging vanishes. However, strange 
stars may possess independent electric dipole and neutron stars independent electric quadrupole 
moments that may lead to the presence of purely electromagnetic contributions to the frame dragging. 
Moreover, recent observations have shown that in stars with strong electromagnetic fields, the 
magnetic quadrupole may have a significant contribution to the dynamics of stellar processes. As an 
attempt to characterize and quantify the effect of electromagnetic frame-dragging in these kind of 
astrophysical objects, an analytic solution to the Einstein-Maxwell equations is constructed here on 
the basis that the electromagnetic field is generated by the combination of arbitrary magnetic and 
electric dipoles plus arbitrary magnetic and electric quadrupole moments. The effect of each multipole 
contribution on the vorticity scalar and the Poynting vector is described in detail. Corrections on 
important quantities such the innermost stable circular orbit (ISCO) and the epicyclic frequencies 
are also considered. 


I. INTRODUCTION 

Frame dragging (Lense-Thirring effect) is the 
quintessential hallmark of general relativity and is the 
result of the spacetime vorticity. It was detected by the 
Gravity Probe B [1] and traditionally, it has been asso¬ 
ciated with rotating stellar astrophysical objects and in 
other astrophysical contexts, such as Galactic models, the 
frame dragging has been associated to the presence of 
magnetogravitational monopoles [2, 3]. Surprisingly, if 
the astrophysical object does not rotate but possesses 
both electric and magnetic fields, the spacetime vorticity 
does not vanish [4-6]. In this case, frame dragging is of 
purely electromagnetic nature and it is associated with 
the existence of a nonvanishing electromagnetic Poynting 
vector around the source [5-9]. 

In most of the early models of astrophysical objects, 
macroscopic charge neutrality is assumed [10] and the 
magnetic field is characterized in terms of a pure dipole 
component [11-14]. This idea was endorsed by the con¬ 
firmation of many predicted features of a star with a 
dipole field using two- and three-dimensional magneto- 
hydrodynamic numerical simulations of magnetospheric 
accretion [15-18]. Thus, under this configuration the 
purely electromagnetic contribution to the vorticity ten¬ 
sor vanishes. However, astrophysical objects such as 
strange stars may possess independent electric dipole (see 
Sec. 10.4 in Ref. [19]) and neutron stars independent elec¬ 
tric quadrupole moments (see below), which may lead 
to the presence of purely electromagnetic contributions 
to frame dragging. A precise account of those contribu¬ 
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tions is not available because of the lack of an analytic 
exact solution with such a complex electromagnetic field 
configuration. 

Moreover, recent observations have shown that in 
stars with strong electromagnetic fields, the magnetic 
quadrupole may have significant contributions on the dy¬ 
namics of stellar processes. Therefore, it is well known 
by now that the actual configuration of the magnetic 
field of strongly magnetized stars may depart from the 
dipole configuration [20]. As an attempt to describe the 
spacetime geometry surrounding these kind of astrophysi¬ 
cal objects, an analytic solution to the Einstein-Maxwell 
field equations is constructed here on the basis that the 
electromagnetic field is generated by the combination 
of arbitrary magnetic and electric dipole plus arbitrary 
magnetic and electric quadrupole moments. This ana¬ 
lytic exact solution allows for analyzing, e.g., the effect 
of each multipole contribution on the vorticity scalar and 
the Poynting vector (see Sec. V). Moreover, it is possi¬ 
ble to predict corrections to important quantities such 
as the innermost stable circular orbit and the epicyclic 
frequencies. 

To motivate further the derivation of the model consid¬ 
ered here, the observational evidence for the existence of 
nondipolar fields in a variety of astrophysical objects is 
discussed next. 


II. OBSERVATIONAL EVIDENCE OF 
NON-DIPOLAR FIELDS 

Measurements of magnetic fields of strongly magnetized 
stars, based on the Zeeman-Doppler imaging technique 
[21], have shown that for these kind of astrophysical ob¬ 
jects the magnetic field has a complicated multipolar 
topology in the vicinity of the star [22-26] . This feature 
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certainly is of prime relevance in, e.g., the accretion-disk 
dynamics in binary systems because if the quadrupole 
component dominates, then the flow of matter into the 
star will certainly differ from the well-known dynamics 
induced by a pure dipole field [25]. 

Complex configurations of magnetic fields are also 
present in stars such as the T-Tauri stars. They are 
young stellar objects of low mass that present variations 
in their luminosity. An important subclass of this kind 
of stars are the so-called classical T-Tauri stars (cTTS) 
because they present accretion from the circumstellar disk 
[27]. Recent evidence points out that in classical T-Tauri 
stars, the magnetic field near the star is strongly non- 
dipolar [26]. In the particular case of V2129 Oph, there 
exits a dominant octupole, 0.12 T, and a weak dipole 
component, 0.035 T, of the magnetic field [28]. Under¬ 
standing the circumstellar disk dynamics, under complex 
field topologies, could provide insight into the formation 
of planets and the evolution of the star itself. 

The discussion above also applies to white dwarfs. The 
first observations indicated that only a small fraction of 
white dwarfs appear to exhibit magnetic fields. However, 
the observational situation changed significantly by the 
discovery of strong-field magnetic white dwarfs [29-31], 
which are known to cover a wide range of field strengths 
~ 1 — 100T and deviates from the simple dipolar configu¬ 
ration [32, 33]. Recent spectropolarimetric observations 
in white dwarf have shown that, in addition to the dipole 
term, the quadrupole and octupole terms make significant 
contributions to the field when it is represented as an 
axisymmetric multipolar expansion [34]. Many studied 
cases indicate that higher multipole components or non 
axisymmetric components may be required in a realistic 
model of white dwarfs (see, e.g., Ref [35]). 

Magnetars constitute an additional source of motiva¬ 
tion. They are characterized by their extremely powerful 
magnetic fields, covering strengths from ~ 10 s to ~ 10 11 T 
[36, 37], so that they can have occasional violent bursts. 
However, certain magnestars, such as the SGR 0418+5729, 
undergo this bursting phenomenon even with a weak mag¬ 
netic fields (~ 10 8 T) [38], For SGR 0418+5729, the 
observed X-ray spectra cannot be fit with this low field 
strength; hence, it has been suggested that a hidden non¬ 
dipole field component must be present to explain the 
bursting episodes [39]. 

In summary, there is sufficient observational evidence 
to develop a consistent analytic closed representation 
of the exterior spacetime around stars with non-dipolar 
magnetic fields. 


III. ANALYTIC FORMULAE OF THE MODEL 

By combining the facts that (i) almost all the analytic 
closed form models for relevant astrophysical objects have 
been conceived in the frame of stationary axisymmetry 
geometry (see [40-45] for the case of neutron stars), (ii) 
powerful tools to construct exact solution to the Einstein- 


Maxwell field equation have been developed, e.g. [46-48], 
and (iii) systematic studies on the construction of exact 
solution from its physical content have been performed 
[49, 50], a new analytic exact solution to the Einstein- 
Maxwell field equations is introduced below. This so¬ 
lution provides physical insight, e.g., into the influence 
of high order electromagnetic multipole moments in the 
frame dragging and in studying quasi-periodic oscillations 
(QPOs), which become a useful tool to identify the char¬ 
acteristics of the compact objects present in Low Mass 
X-Ray Binaries (LXRB) [51-53]. The model presented 
here is a member of the iV-solitonic solution derived in 
Ref. [54]. In Appendix A the relevant equations of the 
derived metric are summarized. 

In terms of the quasicylindrical Weyl-Lewis-Papapetrou 
coordinates x M = (t,p,z,<j>), the simplest form of the line 
element for the stationary axisymmetric case was given 
by Papapetrou [55], 

ds 2 = g^dx^dx", (1) 

with g tt = -f(p,z), g t c/> = f{p 1 z)uj{p,z), g^ = 
P 2 f~ 1 {p,z) - f(p,z)uj 2 (p,z) and g zz = g pp = 
e 2 7(p>U/ _1 (p, z). The metric functions /, to and 7 can 
be obtained from the Ernst complex potentials f and $ 
(see details in Ref. [46]). The Ernst potentials obey the 
relations [46] 

(Ref + |<L| 2 )V 2 f = (Vf + 2$*V4>) • Vf, 

(Ref + |$| 2 )V 2 $ = (Vf + 24>*V<L) • V$ . 

From a physical viewpoint, the Ernst potentials are rele¬ 
vant because they lead to the definition of the analogues of 
the Newtonian gravitational potential, £ = (1 — f )/(l + f), 
and the Coulomb potential, q = 24>/(l +f). The real 
part of £ accounts for the matter distribution and its 
imaginary part for the mass currents. Besides, the real 
part of the q potential denotes the electric field and its 
imaginary part the magnetic field. 

The Ernst equations (2) can be solved by means of the 
Sibgatullin’s integral method [47, 48], according to which 
the complex potentials f and 4> can be calculated from 
specified axis data f (z,p = 0) and &(z,p = 0) [47, 48]. 
Motivated by the accuracy [41, 42] and the level of gen¬ 
erality of the analytic solution derived in Ref. [42], the 
Ernst potential £{z,p = 0) is chosen as in Ref. [42], To 
construct the exact solution that represents the electro¬ 
magnetic field configuration described above, the Ernst 
potential 4>(z,p = 0) is chosen following the prescription 
in Ref. [49]. The Ernst potentials on the symmetry axis 
read 

z 3 - z 2 (m + ia) — kz + is 

£(z,p = 0) = 3 2 , -rx— . , 

z A + z z [m — ia) — kz + is 

»(*,» = 0 )= Z + H V + i Vk + + X ■ 

z A + z z (m — m) — kz + is 

The physical meaning of the parameters in Eq. (3) is 
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derived from the multipole moments calculated using 


the Fodor-Hoenselaers-Perjes procedure [56] (see also 
Ref. [57]). For the present case, 


Pq = to, P\ — iam , P 2 = (k — a 2 )m, P 3 = —im(a 3 — 2 ak + s ), 

P 4 = [70a 4 — 210a 2 k + 13 a<;(n — iv) + 140as + 10fc(7fc — m 2 + ^ 2 ) + 3(/z 2 + v 2 )] , 

P 5 = — 21 21a 5 + 84 a 3 k + a 2 (— 6 /x? — 63s + 6 ivQ 

+ a [—63fc 2 + 6 k(m — ?)(m + ?) + fj 2 + 5+ v 2 + 5%<r] 

+ k (3/z? + 42s - ivQ) + 2i(n( + v\) + 7s (? 2 - to 2 )} , 

Qo = <T, Qi = V + *(ac + /x), Q 2 =—a 2 <; - an +k<; + x +i(av + (), 

Qo = —a 2 v — aC, + kv + i(— a 3 ? — a 2 n + a(2k<; + x) + k/i — sc), 

Q 4 = a 4 c + a 3 (/x — iv) + a 2 (—3fcc — iC, — x) + [l40sc — (n — u>)(140fc — 3m 2 — 10 c 2 )] 
+ y {c [c(fcc + + x) + (m - A;) 2 ] + (7 k - m 2 )(ki; + i( + x) + 7s(/x - iv)} , 

Q 5 = ^- |21ia 5 c + 21a 4 (u + in) — 21m 3 (4fc? + + x) + ° 2 [63isc — i(63/c + 8 c 2 )(n — w)] 

+ ia {—c [c(B/cc + + X) + 9 ^ 2 — 16*/zu — 7u 2 ] + (21/c — 2 to 2 )(3/cc + 2i£ + 2%) 

+ 42s(^i — iv)} — i(n — iv) (— 21 k 2 + 2 km 2 + ^i 2 + v 2 ) + k<; 2 (v — in) 

+ c(—42 iks — 6 nC + 6 ?’/xx + 7 im 2 s + 8 iC,v + 8 vx) + 21s(£ — ix ) + 7*sc 3 } • 


is the usual rotation-induced deformation and a second 
contribution km that accounts for a possible intrinsic de¬ 
formation of the star [59, 60]. An analogous argument can 
be formulated in the case of the electric moments. The 
real part of Q - 2 accounts for the electric quadrupole contri¬ 
bution to the total electromagnetic quadrupole moment. 
The terms —a 2 c and —an account for the rotation-induced 
redistribution of the electric charge and deformation of 
the magnetic dipole; whereas, the term fcc accounts for 
the contribution from the charge distributed over the 
intrinsic deformed mass. The additional parameter \ is 
added to account for any additional possible contribution 
to the total quadrupole moment. 

The multipole expansion in Eq. (5) shows that even if 
the magnetic dipole parameter is zero (n = 0 ), a magnetic 
dipole component (imaginary part of Qi) is present pro¬ 
vided by the rotation of the electric charge Qo- Similarly, 
even if the magnetic quadropole parameter £ is set to 
zero, a rotating electric dipole can induce a magnetic 
quadrupole (imaginary part of Q 2 ). For the electric part 
of the multipole expansion in Eq. (5), an analogous be¬ 
havior is observed, namely, a rotating magnetic dipole 
can induce an electric quadropole moment and a rotating 
magnetic quadrupole can generate an electric octupole 
moment. Based on these processes, the astrophysical 
source can afford a non-vanishing induced Poynting vec¬ 
tor and, correspondingly, an induced non-vanishing flux 
of electromagnetic energy around the source that will 


Symbol 

Associated Multipole Moment 

? 

Electric monopole 

V 

Electric dipole 

X 

Electric quadrupole 


Magnetic dipole 

C 

Magnetic quadrupole 


Table I. Summary of the electromagnetic parameters and the 
multipole moments they are related to. 

Specifically, the interpretation of the parameters based on 
the multipole expansion in Eqs. (4) and (5) is as follows. 
The real parameter m corresponds to the total mass, a 
to the total angular moment per unit mass while k and 
s are related to the mass-quadrupole moment and the 
differential rotations, respectively. For later convenience, 
an electric monopole contribution Q 0 , characterized by 
the parameter q, was introduced above. Parameters v 
and n are associated with the electric and magnetic dipole 
moments, respectively, whereas x and £ with the electric 
and magnetic quadrupole moments, respectively. The 
existence of an electric dipole is theorized for strange stars 
(See Ref. [58]). A summary of the arbitrary parameters 
and the multipole moments they are related to can be 
found in Table I. 

The mass moment P 2 governs the deformation of the 
star and it is composed of two parts: the term a 2 m that 
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contribute to the frame dragging induced solely by the 
mass currents. 


IV. CHARACTERIZATION OF THE 
ELECTROMAGNETIC FIELDS 


To describe the electromagnetic properties of the so¬ 
lution, the electric and magnetic fields, in the spacetime 
surrounding the star, are calculated by means of the 
expressions 

E a = F a/3 uP, B a = ~ t a p' 1& F 1 su p , ( 6 ) 

where F a p is the electromagnetic field tensor F a p = 
2 A [ 0 ;a ], = ( 0 , 0 , A^,— A t ) is the electromagnetic four- 

potential, u a is a time-like vector and e Q / 3 7< 5 is the totally 
antisymmetric tensor of positive orientation with norm 
= —24 [61]. For a congruence of observers at 
rest in the frame of ( 1 ), the four-velocity is defined by 
the time-like vector u a = (1/y/J, 0,0, 0). The vectorial 
fields have components in the p and z directions only. 
The components of the electric field are given by 


E P = 


77 

e 2 7 


At,p ) 


E z 


77 

e 2 7 


A 


t,Z 5 


( 7 ) 
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Figure 1. Magnetic field force lines for m = 2.071 km, 
j = 0.194, Q = -2.76 km 3 , s = -2.28 km 4 , ? = 0 km, 
v = 0 km 2 , x = 0, p, — 1 km 2 for (a) £ = 0 km 3 , (b) 
£ «= 15 km 3 , (c) £ = 30 km 3 and (d) £ = 50 km 3 . The non- 
electromagnetic parameters correspond to the model 2 for the 
equation of state L in Ref. [62] (see also Table II). 


and for the magnetic field by 

f 3/2 

B p = — 2^7 (— 1 ^*,2 + A ^ z ), ( 8 ) 

pe^ 

f3/2 

Bz = + A(j>,p) ■ (9) 

pe 

The explicit form of the fields can be found in Appendix B. 


Model 

to [km] 

j 

Q [km 3 

s [km 4 

M2 

2.071 

0.194 

-2.76 

-2.28 

M3 

2.075 

0.324 

-7.55 

-10.5 

M4 

2.080 

0.417 

- 12.2 

- 22.0 

M5 

2.083 

0.483 

-16.2 

-33.9 


Figure 1 shows the force lines of the magnetic field for 
various values of the magnetic quadrupole parameter £ 
and for realistic values of the mass and mass current 
multipoles. Specifically, the vacuum multipole moments of 
the solution mass, angular moment and mass quadrupole 
and current octupole have been fixed to the numerical 
ones obtained in Ref. [62]. They are listed in Table II. 
In particular, for Fig. l.a, £ = 0 km 3 ; l.b, £ = 10 km 3 ; 
l.c, £ = 25 km 3 and for Fig. l.d, £ = 50 km 3 . The 
increasing of the separation between consecutive force 
lines indicates that the magnetic field decreases while 
the distance increases. Figure 1 not only shows how the 
reflection symmetry around the plane z = 0 is broken 
because of the magnetic quadrupole [ 6 ], but also shows 
that at large distances from the source, the magnetic field 
behaves like a magnetic dipole despite the presence of 
strong non-dipolar contributions. 

Figure 2 shows the force lines of the electric held for a 
variety of values of the electric quadrupole for realistic 
values of the mass and mass current multipoles listed 


Table II. Realistic numerical solutions for rotating neutron 
stars derived by Pappas and Apostolatos [62] for the equation 
of state L. Here, m is the total mass of the star, j is the 
dimensionless spin parameter: j = J/M 2 (being J the angular 
momentum), Q is the quadrupole moment and s is the current 
octupole moment. See Table VI of [62]. 


in Table II. Specifically, for Fig. 2. a, y = 0 km 3 ; 2.b, 
X = 10 km 3 ; 2.c, y = 25 km 3 and for Fig. 2.d, y = 50 km 3 . 
As above, it is clear that reflection symmetry is broken 
and that at large distances the electric dipole component 
dominates the held conhguration. 

In general, the electromagnetic held of astrophysical 
objects is expected to be a combination of the results 
depicted in Figs. 1 and 2. Moreover, as shown below, 
the breaking of the rehection symmetry has an important 
role on the vorticity of the spacetime. After character¬ 
ize the electromagnetic helds, the generation of purely 
electromagnetic frame dragging is discussed below. 
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Figure 2. Electric field force lines for m = 2.071 km, j = 0.194, 
Q = —2.76 km 3 , s = —2.28 km 4 , ? = 0 km, v = 1 km 2 , 
p = 0 km 2 , ^ = 0 for (a) x = 0 km 3 , (b) x = 10 km 3 , (c) 
X = 25 km 3 and (d) x = 50 kin 3 . 


as discussed below, it could decrease the total vorticity 
of the spacetime. Figure (3) depicts the Poynting vector 
circulation around the source, when /i is chosen parallel 
(l.h.s. panel) and antiparallel (r.li.s. panel) to the star’s 
rotation axis. 



-40 -20 0 20 x[km] -40 -20 0 20 i[km] 


Figure 3. Vector field of the Poynting vector for the star model 
m = 2.071 km, j = 0.194, Q = -2.76 km 3 , s = -2.28 km 4 . 
The electromagnetic multipoles were chosen zero except c = 
0.1 km, p = 1 km 2 in the left panel and p = —1 kin 2 in the 
right side. The direction of the Poynting vector is due the 
sign change in the multipole expansions in Eq. (5). Figures 
are pr esented in the quasi-Cartesian auxiliary coordinates 
r —> \/x 1 + y 2 and <j> = arctan (y/x). 
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V. VORTICITY SCALAR AND POYNTING 
VECTOR 

The physics of the frame dragging states that the rota¬ 
tion of the source induces a twist in the neighbourhood 
that drags any frame of reference near the source. In the 
case of spacetimes endowed by complex electromagnetic 
fields and mass currents, frame dragging originates from 
a combination of the vorticity of the electromagnetic field 
and the vorticity associated with the mass currents. 

Poynting vector .-The Poynting vector S carries the 
information about the electromagnetic energy flux in the 
spacetime. Because of the axially symmetric character 
of the spacetime, only the component along the unitary 
vector e^ survives [6]. In terms of the Ernst potentials 
S^ = 7/Im /(47rpe 27 ) [63] or more conveniently 

S 0 = ^ V$* x V<&, (10) 

8t rpe 27 v ' 

where it is clear that S& vanishes if $ is purely real, purely 
imaginary or when their real and imaginary parts are pro¬ 
portional to each other. Due to the complex combination 
of the electromagnetic moments with the mass currents, 
their effects on the Poynting vector (and subsequently to 
the vorticity) are often subtle. However, the multipole ex¬ 
pansion in Eq. (5) allows for a detailed description of each 
contribution. For instance, if v = 0, x = 0 and |ac| < |/i|, 
it is then clear that, to leading order, a sign change in the 
magnetic dipole moment parameter p changes the sign 
of the magnetic field. This has the effect of changing the 
rotation direction of the Poynting vector (see Fig. 3) and 


The change of the circulation of the Poynting vector can 
be undesrtood from the potential q(p, z) that is analoguos 
to the Coloumb potential (see above). Under the condition 
that the reflection symmetry is not broken, v = 0 and 
£ = 0, a sign change of the magnetic dipole parameter /i 
is equivalent to change the sign of the coordiante z with 
a global minus sign, namely, q(p, z: v = 0, C = 0, —a) = 
= 0, C = 0,p) = -Hp^v = 0, C = o’p). 
Therefore, under the conditions that reflection symmetry 
imposses [49], all the electric field moments change sing 
and this leads to a global sing of the Poynting vector. 

Vorticity Scalar- For an observer at rest respect to (1), 
the vorticity tensor is defined by u> a p = it[ a; / 3 ] + U[ a upy 
The direct calculation of u a p for (1) yields to [6] 


/ 0 0 0 

0 0 0 

0 0 0 

\o Wf^, P Wfu, z 


0 \ 

-hVfu,p 

VJv,z 

0 / 


( 11 ) 


To quantify the vorticity, it is convenient to introduce the 
vorticity scalar that is defined by the contraction of the 
vorticity tensor and reads 

w v = (ufgwf) 2 = fyJf{u* p + uP z )ly/2e 2 ~ip 2 , (12) 


with = — pf 2 3(f z + 2$*$ i2 ) and w 2 = pf 2 3(£ p + 
24>*$ p). More explicitly, it reads 


U v = 


o—7 


72 / 




2<F*$ 




2<I>*$ p ] 2 


(13) 
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The vorticity scalar can be understood in terms of its 
fluid mechanics analogue. It represents the rotation of 
the fluid. In the general relativistic case, it can be related 
to the rotation velocity of a family of congruences. 

Equation (13) is very useful to analyze the contribu¬ 
tions to the vorticity because the electromagnetic and 
mass currents terms can be easily identified there. More¬ 
over, this identification can be accompanied by further 
expressing these contributions in terms of the parameters 
of the Ernst potentials (3). In particular, in absence of 
electromagnetic fields, the imaginary part of £ is asso¬ 
ciated with mass currents. Thus, for static sources one 
could attempt to assume £ as real in Eq. (13) and focus 
only in the electromagnetic contributions encoded in $. 
However, because £ and <f> are not independent objects 
[see, e.g., Eq. (2)], in the presence of electromagnetic fields, 
£ and $ are complex even for non-spinning astrophysical 
objects. Notwithstanding, if <f> is purely real or purely 
imaginary, for non-rotating objects, the imaginary part 
of £ vanishes [see, e.g., Eq. (2)]. Therefore, interest here 
is in identifying when the electromagnetic contribution 
does not vanish. 

In doing so, it is convenient to express w v = |3(V£ + 
2<h*V$)|/(v / 2Je 7 ), so that the term vanishes if 

$ is purely real, purely imaginary or when their real and 
imaginary parts are proportional to each other. Based 
on the discussion above, see also Eq. (5), the purely 
real case corresponds to the absence of magnetic fields 
and the purely imaginary case to the absence of electric 
fields. In the cases when the Poynting vector is zero, the 
electromagnetic contribution to the vorticity vanishes; this 
generalizes the results in Ref. [6] where a particular space- 
time was considered. The particular case of proportional 
real and imaginary parts leads to the case of proportional 
electric and magnetic fields and corresponds to the cases 
studied in Refs. [4, 6]. 

Note that the considerations above on the vorticity 
scalar and the Poynting vector are completely general 
and valid for any stationary axially symmetric spacetime. 

To study the contribution of the electromagnetic energy 
flux to the vorticity, characterize first the contributions of 
the rotation of the source. In this case, the Ernst potential 
<!> that encodes all the electromagnetic moments is zero 
and £ is complex. For fast rotating stars, the vorticity 
scalar (and the Lense-Thirring effect) is larger than for 
slow rotating stars, this can be seen in the Fig. (4). 

When one includes an electromagnetic field to a partic¬ 
ular star model, the Ernst potential <1> is no longer zero, 
and it has an important role in the vorticity scalar. The 
particular contribution depends on the structure gener¬ 
ated by the multipole expansion and the Poynting vector. 
Before considering the purely electromagnetic contribu¬ 
tion to the frame-dragging in realistic situations (zero 
or negligible total electric charge), consider the case of a 
source with a dipole magnetic field and electric charge. 

The upper panel of Fig. (5) depicts the functional de¬ 
pendence of the vorticity scalar on the distance from the 



Figure 4. Vorticity generated by the mass currents. All the 
electromagnetic moments are zero in this case. The vorticity 
of the spacetime was calculated for the models 2 (j = 0.194), 3 
(j = 0.324), 4 (j = 0.417) and 5 ( j = 0.483) with the equation 
of state L in Ref. [62]. They are listed in the Table II 

star for a fixed value of the electric monopole and for a 
variety of magnetic dipole moments. The vorticity scalar 
decreases (increases) from its value in the vacuum situa¬ 
tion (p = 0) when the dipole is parallel (anti-parallel) to 
the star’s rotation-axis. The reason for this phenomenon 
can be understood in terms of Eq. (13). To do so, note 
that £ decreases monotonically as a function of ^ and p\ 
thus, its derivatives carries a negative sign. Although the 
same argument applies to 4>, when the sign of p changes, 
see above, the Poynting vector changes its circulation 
and that is enough to change the sign of its contribu¬ 
tion to the vorticity scalar [see Eq. (13)]. Alternatively, 
assume that the reflection symmetry exits and set the 
mass currents parameters a and s to zero, for this case, 
&(p,z-,a = 0, s = 0 ,—p) = 4>*(p, z;a = Q,s = 0,p) and 
therefore V4? changes sign. 

The lower panel of Fig. 5 depicts the functional depen¬ 
dence of the vorticity scalar on the distance from the star 
for a fixed value of the electric dipole and for a variety 
of magnetic dipole moments. The contribution in this 
case is weaker than in the case of an electric monopole. 
However, results in the lower panel of Fig. 5 are more 
realistic than those in the upper panel. Interestingly, in 
the lower panel, there is no change in the sign of the 
electromagnetic contribution to the vorticity scalar when 
p changes sign. The reason for this relies on the fact that 
for v ^ 0 or £ ^ 0, the reflection symmetry around the 
equatorial plane breaks down and the arguments provided 
above do not apply, i.e., <1?(p, z; a = 0, s = 0, —p) does not 
equate to 4>*(p, z; a = 0, s = 0, p). 

Based on the processes described above on the genera¬ 
tion of a Poynting vector from the fields induced by the 
mass currents (a / 0 and s ^ 0), it is clear that even 
if the parameters associated with the electric field (?, v 
and x) are set t° zero, but the parameters associated 
with the magnetic field (p and £) are non-vanishing, then 
a rotationally-induced electromagnetic contribution to 
the vorticity scalar is present. A similar scenario takes 
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Figure 5. Electromagnetically generated vorticity by the pres¬ 
ence of a magnetic dipole and a fixed electric monopole of 
? = 0.1 km (upper panel with w = 0, y = 0 and £ = 0) and an 
electric dipole of w = 10 km 2 (lower panel with ? = 0, x = 0 
and £ = 0). The vorticity of the spacetime was calculated for 
the model 2 for the equation of state L in Ref. [62]. Parameters 
are listed in Table II. 

place if the parameters associated with the magnetic field 
are set to zero and non-vanishing electric parameters are 
considered. In particular, even in the low rotation regime, 
a non-negligible electric field is generated [64] and may be 
important to characterize the evolution of the electromag¬ 
netic structure in neutron stars. Moreover, in the case of 
fast rotation, the frame-dragging caused by a Kerr black 
hole significantly distorts the structure of an external 
magnetic field and this scenario may be relevant for low 
accreting black holes as the one present in the Milky Way 
center [65, 66]. 

The most relevant contributions to the vorticity scalar 
from the electromagnetic field were considered above. 
However, in general, all the multipole moments contribute 
to the vorticity scalar and the details of the net result may 
deviate from those discussed above; albeit, the magnitude 
of the effect is not expected to differ from the predictions 
in Fig. 5. 


VI. ORBITAL EQUATORIAL MOTION AND 
EPICYCLIC FREQUENCIES 

In general relativity, as in Newtonian gravitation, all 
the physical characteristics of the source have an effect 


on the dynamics of orbiting of particles. One of the most 
distinctive points between both theories is the frame¬ 
dragging, but is not an observable by itself. A way to 
measure it, is to appeal to the dynamics around the 
source and characterizing the behavior of, e.g., neutral 
test particles [40, 53]. 

Observationally, the Keplerian motion of matter could 
be useful to model the quasi-periodic oscillations (QPOs) 
that are present in the low-mass X-ray binaries (LMXBs) 
containing a neutron star [67]. These oscillations occur 
at frequencies in the range of kHz and come in pairs, the 
upper and the lower mode correspond to the frequencies 
of Keplerian motion and periastron precession of the 
accreted matter in the close vicinity of the star [52], An 
additional effect is that these equatorial orbits will exhibit 
a relativistic nodal precession due to frame dragging [53], 
causing a detectable signal in the spectra of LMXBs [68] . 

The following subsections are devoted to quantifying the 
magnetic quadrupole effect in the ISCO and the epicyclic 
frequencies for different models of neutron stars [62]. 

Influence of the field on the ISCO radii 

The dynamics of an orbiting particle can be analyzed 
by using the Lagrangian formalism. To do so, consider 
a particle of rest mass mo = 1 moving in the space- 
time described by the metric functions in Eq. (1), the 
Lagrangian of the particle is then given by 

£ = (14) 

where the dot denotes differentiation with respect to the 
proper time r, :e m (t) are the Weyl-Lewis-Papapetrou 
coordinates. For a stationary and axisymmetric spacetime, 
there are two constants of motion related to the time 
coordinate t and azimuthal coordinate <fi . Thus, the energy 
and the canonical angular momentum, respectively, are 
conserved (see Ref. [69] for details). Assuming that the 
motion takes place in the equatorial plane of the star z = 0, 
and taking into account the four velocity normalization 
for massive particles g^u^u 11 = — 1, one can identify an 
effective potential that characterizes the motion in the 
plane (see, e.g., Ref. [70]) 

V eS (p) = l- E29 ** + 2 2EL9t * + L29tt . (15) 
9uf> 94‘4>9tt 

For a particle moving in a circular orbit, the energy E 
and the canonical angular momentum L are determined 
by the conditions U e ff(p) = 0 and dI4ff/d p = 0 (see, 
e.g., Ref. [71]). The marginal stability condition reads 
d 2 Kff/dp 2 = 0. Thus, the ISCO’s radius is determined 
by solving numerically the previous equation for p. 

The importance of determining the ISCO is that ac¬ 
cretion disks extend from the last stable orbit to exterior 
zones, so the ISCO is an inner boundary for the accreted 
matter. As stated above, the inclusion of a magnetic 
























field causes a dependence of all physical quantities of the 
intensity on the field. From the International System of 
Units (SI), the conversion to geometrized units is given 
by 


A^geom 


10 6 \/Gfi o 
- 2 -MSI) 


(16) 


where G is the gravitational constant, go is the vacuum 
permittivity, c is the speed of light and p is given in units 
of Am 2 . The units of £ are Am 3 and the conversion for 
the quadrupole reads 


Cgeom — 


10~ 9 VWo 


Csi- 


(17) 


In the same way, the electric multipole moments can be 
written in geometrized units as: 


Sgeom — 



^geom 


Xgeom 


10” 3 

nr 6 



Gp 0 

477 


ASI- 


(18) 

(19) 

( 20 ) 


Consider a superposition of a fixed dipolar and 
quadrupole component that will be varied from 1 km 3 to 
50 km 3 . Figure 6 depicts the ISCO radius as a function of 
the parameter p and £ for three realistic numerical solu¬ 
tions for rotating neutron stars models derived in Ref. [62] . 
The parameter p is set to 1 knr 2 that corresponds to a 
magnetic dipole field of 10 12 T. The quadrupole parame¬ 
ter £ is chosen between 0 and 50 km 3 that corresponds 
to magnetic quadrupoles from 0 to 5 x 10 35 Am 3 , respec¬ 
tively. In all these cases, the ISCO radius decreases for 



Figure 6. ISCO radius as a function of magnetic quadrupole 
parameter £, with p = 0 km 2 . The physical parameters for 
the star correspond to the models 2-4 listed in Table II. We 
can se that the ISCO radius decreases for increasing £. 

increasing £). This can be explained as a result of the 
deformation of the space-time by the energy stored in the 
electromagnetic (see, e.g., Ref. [40]). We elaborate more 


on this in the next section. 


Influence of the electromagnetic field in the 
Keplerian and the epicyclic frequencies 


Imposing the conditions of constant orbital radius, 
dp /dr = 0 and taking into account that dip/dr = 
Oicdf/dT, the Keplerian frequency reads (see, e.g., 
Ref. [69]) 


d<p gtc/),p ± \J{gt<p, P )~ 

“K = -T7 = -j (41) 

dt ff00,p 


where “+” and ” denotes the Keplerian frequency for 
corotating and counter-rotating orbits, respectively. In 
Ref. [40] the functional dependence of the Keplerian fre¬ 
quency on the magnetic dipole parameter p was discussed, 
here interest is in the functional dependence on the mag¬ 
netic quadupole parameter The upper panel of Fig. 7 
depicts the functional dependence on £, as in the case 
analyzed in Ref. [40], the Keplerian frequency increases 
with increasing £ because the ISCO’s radius decreases 
(see Fig. 6) as a consequence of the additional deforma¬ 
tion of the space-time by the electromagnetic field. In 
particular, the energy of the electromagnetic field shifts 
the marginally unstable region toward the star. 

Analytical expressions for the radial and vertical fre¬ 
quencies follow from allowing slightly radial and vertical 
perturbations of the orbit. According to Ref. [69], the 
radial and vertical epicyclic frequencies are given by 


27T 


(gu + gtM 2 [ 


2(gtt + 9t^k){gt4 + gu^k) ( 2 


VP" 


( 22 ) 


+(gt 4 + g^^k ) 2 


with a = {p, z}. In the relativistic precession model 
(RPM) [51], the periastron jT and the nodal frequencies 
are the observational relevant ones, they are defined by 


/P = 


Uk 

2-7T 


(23) 


At the ISCO, the radial oscillation frequency equals the 
Keplerian frequency and only the vertical precession 
a = z is considered here. The influence of the mag¬ 
netic quadrupole on the vertical precession frequency is 
depicted in the lower panel of Fig. 7. The behaviour and 
the underlying physical mechanism are analogous to those 
for the Keplerian frequency. 

The frequency i/lt that characterizes the Lense- 
Thirring effect is given by (see, e.g., Ref. [69]) 


i/ LT = — 


1 Stcfi 
g^ 


(24) 
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Figure 7. Influence of the magnetic quadrupole moment £ in 
the epicyclic frequencies for different models of neutron stars 
with ? = 0, [i = 1 km 2 , v = 0 and x = 0. Vacuum parameters 
of the star are listed in Table II. 


and could be relevant to model the horizontal branch oscil¬ 
lations observed in the Low Mass X-Rax Binaries(LMXBs) 
[51]. The strength of the frame-dragging effect increases 
for fast rotating objects and, as shown in the central panel 
of Fig. 7, for stronger magnetic fields. Figure 7 also shows 
that the electromagnetically-induced frame dragging has 
a direct effect on the orbiting particles, even if they are 
neutral. 


VII. CONCLUDING REMARKS 

We present a new stationary axisymmetric nine- 
parameter closed-form analytic solution that generalizes 
the Kerr solution with arbitrary mass-quadrupole mo¬ 
ment, octupole current moment, electric and magnetic 
dipole and electric and magnetic quadrupole moments. 
The analytic form of its multipolar structure and their 
electric and magnetic fields are also presented. Accord¬ 
ing to the arguments presented through the paper, this 
solution could be used to model the exterior gravitational 
and electromagnetic fields around strongly magnetized 
stars, in particular white dwarfs. Besides, this model also 
could be used even for the description of exotic stars such 
as the r Sco recently reported by Donati et al. [72] . 

This solution allowed for a comprehensive analysis 
of the contribution of complex and intense electromag¬ 
netic fields to quintessential hallmark of general relativity, 
namely, the Lense-Thirring effect. It was shown that if 
the value of the parameters is such that the reflection 
symmetry is preserved [49], then a sign change of the 
magnetic dipole parameter fi is enough to change the 
direction of the circulation of the Poynting vector and the 
sign of the electromagnetic contribution to the vorticity 
scalar. The influence of complex electromagnetic fields 
on the ISCO’s radius and the epicyclic frequencies were 
also considered and it was shown that for strong magnetic 
fields, the influence is not negligible. 

Additional interest is in studying the effect of the mag¬ 
netic field in the quadrupole moment of the star (see, e.g., 
Refs. [73, 74]). This topic has gained recently theoretical 
and observational interest and is currently discussed un¬ 
der the title of I — Love — Q relations [74]. Due to the 
great generality of the present analytic exact solution, this 
subject will be explored in a forthcoming contribution 

[75]. 
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Appendix A: Metric Functions of the Developed Solution 


This appendix summarizes the relevant equations of the developed metric. The potentials in the symmetry axis can 
be written as [54]: 


with 
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Then, the Ernst potentials and the metric functions in whole spacetime are derived with the aid of the Sibgatullin’s 
integral method [47, 48]. By using the representation proposed in Ref. [76] and used also Ref. [42], 
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Pii = (-1)*+% ( e l "^23|irTmp + e 2 ^31 |Imnp + £ 3 ^ 121 Imnp) 1 
Pijk = ( — l) 4 +J ’ + 1 Aijfc(eJ ^23|Zmn + 62 ^31|Zmn + e 3 ^12|Zmn) 5 


Ajj; (g^ oij ) Dj , A k (^2 g^) ^fc) -^2 Dj Dk 5 Di . o \/ o \/ 

— pij(di — P2j(di — P 3 J 

= A^2 |mn “1“ -^3(^m) ^12|nZ “I - -^ 3 (^ 22 ) ^12|Zra 5 ^Zra|np = -^z(^n) -^ra(^p) ZZj(dp) ZZ m (d n ) , 

2 ]lp^nK-/ 3 Z*) ftfk 


-^z(^n) 


-*E 


nE(A* - 4*) nLi(A* - a) “ ^ - ajk - a) 

^Zlmnp i?3(az) ^12|mnp ? -^Z|mnp ^Z|mnp “I - ^m|npZ “I” ^nlpZm “I - ^p|Zmr 
^Zm|nps /(^n) ^/m|ps “t - /(^p) ^Zm|sn “t” /(^s) ^Zm|np ; 

^Zmlnprs /(^n) *^/m|prs f (^p) ^/m|rsn d“ f (^r) ^Zm|snp f (^s) ^/m|npr ; 

being a’s the roots of the Sibgatullin’s equation [47, 48] 

e(z) + e(z) + 2 f(z)f(z) = 0 . 


; ^lm\nps ^lm\np 4“ ^Zm|ps 4“ ^lm\s 


(A4) 


Appendix B: Electromagnetic Field: Analytic Form 


The A t potential is the real part of the electromagnetic Ernst potential <f>, and the potential A$ can be calculated as 
the real part of the Kinnersley potential K, = A$ + iA' t [77] , which can be obtained using the Sibgatullin’s method and 
can be written as 


y ■ I (/i 4- /*2) 

K ~-' Trs 


(Bl) 
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Thus, the closed-form expressions for the electric and magnetic fields are 


E P = - 


A 


VW^W+W 

A 


Re<! 


E z = - -=-_ Re<^ 

'/W^W+W I 


C iP -Chi{A-B) if 
A-B 

C tZ -C\n(A- B) >z 
.1 B 


B P = 


Im {A-B)H + (A-B)G-Cl\ A^|A| 2 -|B|2 + |C| 2 j (f 1+ f 2 )[I >z -I]n(A-B) tI 
p\A-B \ 2 z+ p\A-B\ 2 m | A — B 


_ Im(A - B)H + (A - B)G - Cl\ A^\A\ 2 - \B\ 2 + |C| 2 / (/i + / 2 ) [I, p - /ln(A - B), p 

Bz ~ p\A — B\ 2 Kp + p\A~B\ 2 lm ] A^B 


when 


i^i 2 nLi^ 

\A-B\ 


A = 



